Reconciling Semiclassical and Bohmian Mechanics: 
I. Stationary states 



OO 
O 

o 

(N 
m 

Oh' 



> 

m 

O 

OO 

o 



X 



Bill Poirier 

Department of Chemistry and Biochemistry, and Department of Physics, 
Texas Tech University, Box 4IO6I, Lubbock, Texas 79409-106U 

The semiclassical method is characterized by finite forces and smooth, well-behaved trajectories, 
but also by multivalued representational functions that are ill-behaved at caustics. In contrast, 
quantum trajectory methods — based on Bohmian mechanics (quantum hydrodynamics) — are char- 
acterized by divergent forces and erratic trajectories near nodes, but also well-behaved, single- valued 
representational functions. In this paper, we unify these two approaches into a single method that 
captures the best features of both, and in addition, satisfies the correspondence principle. Stationary 
eigenstates in one degree of freedom are the primary focus, but more general applications are also 
anticipated. 



I. INTRODUCTION 

Theoretical and computational physical chemists have 
long sought reliable and accurate means of performing 
quantum dynamics calculations for molecular systems, as 
quantum effects such as tunneling and interference often 
play an important role in such systems. If "exact" meth- 
ods are required — i.e., numerical techniques for which the 
error bars may (in principle) be reduced arbitrarily — 
the traditional approach has been to represent the sys- 
tem Hamiltonian using a finite, direct-product basis set. 
However, this approach suffers from the drawback that 
the scaling of computational effort is necessarily expo- 
nential with system dimensionality. 

Recently, a number of promising new methods have 
emerged that may spell an end to exponential scaling — 
or at the very least, drastically reduce the exponent. The 
latter category includes various basis set optimization 
methodsj^i^i^iii^ which have nearly doubled the number 
of degrees of freedom (DOFs) that may be tackled on 
present-day computers, from about 6 to 10 DOFs. Very 
recently, the first basis-set method to defeat exponential 
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scaling (at least in principle) was introduced^^^ 
method, which uses wavelets in conjunction with a phase 
space truncation scheme, has been applied to model prob- 
lems up to 15 DOFs, and is easily extendible to higher 
dimensionalities — although some technical issues vis-a- 
vis applicability to real molecular systems must still be 
resolved. 

A completely different approach to the exact quan- 
tum dynamics problem may be found in time-dependent 
trajectory methods. Although trajectory methods are 
extremely common in molecular dynamics applications, 
they are almost always classical, quasiclassical, or 
semiclassica l^^'^'^1^^'^^1^^'^'''^^'^^ — i.e., not exact, in the 
sense described above. However, it is possible to perform 
exact quantum dynamical propagation using trajectory- 
based methods. These so-called "quantum trajectory 
methods"^i^i^i2^i^ (QTMs) are based on the hydrody- 
namical picture of quantum mechanics, developed over 
half a century ago by Bohm^i^ and Takabayasi^l using 
even earlier ideas of Madelung-- and van Vleckfi^ 



Trajectory methods of all kinds are appealing, because 
they offer an intuitive, classical-like understanding of the 
underlying dynamics. QTMs are especially appealing, 
however — not only because they ultimately yield exact 
results, but also because they offer a pedagogical under- 
standing of quantum effects such as tunnelingi2£i2i Cu- 
riously, QTMs thus far have not fared so well as semi- 
classical methods, vis-a-vis their treatment of another 
fundamental quantum effect — interference. This issue is 
discussed in more detail below, as it is of central concern 
for the present paper. An in-depth comparison of inter- 
ference phenomena is provided in an intriguing article by 
Zhao and Makri,^- Garashchuk and Rassolo v^^i^° discuss 
the interesting connection between QTM and semiclas- 
sical propagators fi^iiiii^ in the context of Hcrman-Kluk 
initial value representations i^ii^^i^ 

Perhaps the greatest attraction of QTMs, however, 
has been the promise that they may be able to defeat 
exponential scaling. In any event, QTMs have under- 
gone tremendous development since their introduction in 
1999 — most notably within the last year or two. Much of 
the early development centered around accurate evalua- 
tion of spatial derivatives of the wavefunction,32i2ii2i but 
with the introduction of local least-squares fit adaptive 
and unstructured grid techniques r^i^^iSi this difficulty is 
now essentially resolved. This has paved the way for a 
number of interesting applications of QTMs, including 
barrier transmission,?" non-adiabatic dynamics, "^^ and 
mode relaxationi^ Several intriguing phase space gen- 
eralizations have also emerged ;2ZiSi2^i^ii£ of particular 
relevance for dissipative systemsi^'^'^i^ 

On the other hand, QTMs still suffer from one ma- 
jor drawback — they are numerically highly unstable in 
the vicinity of nodes. This "node problem" manifests in 
several different waysi^^i^i (1) infinite forces, giving rise 
to kinky, erratic trajectories; (2) compression/inflation 
of trajectories near wavefunction local extrema/nodes, 
leading to; (3) insufficient sampling for accurate deriva- 
tive evaluations. In the best case, this can result in sub- 
stantially more trajectories and time steps than the cor- 
responding classical calculation; in the worst case, the 
calculation may fail altogether, beyond a certain point in 
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time. 

For many molecular applications (though certainly 
not all), the initial wavepacket is nodeless; however, it 
may develop nodes at some later point in time. More- 
over, from a practical standpoint, nodes per se are not 
the only source of numerical difficulty; in general, any 
large or rapid oscillations in the wavefunction — termed 
"quasinodcs" — can be sufficient to cause the problems 
described above. Such oscillations are very easily formed 
in molecular systems, particularly during barrier reflec- 
tion. Accordingly, several numerical techniques are being 
developed to address the node problem, including the Ar- 
bitrary Lagrangian-Eulerian frame methodj^i^ and the 
artifical viscosity methodi^i^ 

In this paper, we take a different approach to the node 
problem, based on a thorough understanding of the dif- 
ferences and similarities between Bohmian and semiclas- 
sical mechanics. Formally, these two theories share many 
similarities, as was known from the earliest day o^'^''^^ — 
yet in practical terms, semiclassical and quantum trajec- 
tories often behave very differently. For instance, the 
former may cross in position space, but not in phase 
space; the latter do exactly the opposite. For the spe- 
cial case of stationary eigenstates in 1 DOF (the focus 
of the present paper), semiclassical trajectories evolve in 
phase space along the contours of the classical Hamilto- 
nian, whereas quantum trajectories do not move at all. 
For well-behaved potentials, classical trajectories are al- 
ways smooth and well-behaved, but quantum trajectories 
may be kinky and erratic. 

As noted by Zhao and Makri)^ nowhere are the dif- 
ferences between the two methods more profound than 
in the treatment of nodes and interference phenomena — 
which is not even qualitatively similar. In the semiclassi- 
cal treatment, the approximate wavefunction is expressed 
in terms of simple functions that are generally as smooth 
and well-behaved as the potential itself. Oscillations and 
nodes are obtained when different lobes or "sheets" of 
the semiclassical wavefunction come to occupy the same 
region of space — thus giving rise naturally to interfer- 
ence. In contrast, the Bohmian representation of the 
wavefunction is single- valued, and therefore incapable of 
self-interference. Consequently, all of the undesirable — 
even "unphysical" — aspects of quantum trajectories, as 
described in the previous paragraph, are necessary in or- 
der to represent nodes and quasinodes explicitly. 

From a dynamical perspective, the only difference be- 
tween semiclassical and quantum trajectories is the quan- 
tum potential, Q, and it is a remarkable fact that Q 
alone is responsible for all of the very fundamental dif- 
ferences described above. On the other hand, this sit- 
uation is also cause for alarm, for Q turns out to be 
the order correction that is ignored in semiclassi- 
cal treatments — being regarded as "insignificant" in the 
large action limit, in accord with the correspondence 
principle. From the previous discussion, it is clearly 
incorrect to regard Q as insignificant, which seems to 
place the semiclassical approximation — and indeed, the 



correspondence principle itself — in jeopardy. Moreover, 
there are certain unappealing aspects of the semiclassi- 
cal approach — multivaluedness, caustics, phase disconti- 
nuities, etcji^iii^ii^iiiii^ — that simply do not arise in a 
Bohmian treatment. On the other hand, the semiclassi- 
cal approximation is known to be valid in the large ac- 
tion limit — which together with the undesirable features 
of the Bohmian approach as discussed in the previous 
paragraph, seems to call into question the physical cor- 
rectness of the latter. This paradox has been a source of 
concern for some researchers, notably Einstein4^i^i^ 



The primary goal of the present paper is to reconcile 
the semiclassical and Bohmian theories, in a manner that 
preserves the best features of both, and also satisfies the 
correspondence principle. At least within the context 
of stationary eigenstates in 1 DOF, the above paradox 
turns out to be remarkably easy to resolve. It can be 
shown that the disturbing dissimilarities described above 
stem not from the theoretical methodologies themselves, 
but from the simple fact that each method uses a dif- 
ferent ansatz to represent the wavefunction — thus giving 
rise to a rather unfair comparison. In particular, since 
the semiclassical functions are double- valued in the clas- 
sically allowed region of space, the stationary wavefunc- 
tion is represented as the sum of two terms — essentially a 
pair of "traveling waves," headed in opposite directions. 
In contrast, the standard Bohmian approach uses a single 
term to represent the wavefunction. Virtually all of the 
disparities described above arise from this simple fact. 



It is therefore natural to consider what would happen 
if the Bohmian formalism were applied to a two-term 
wavefunction — thus placing it on a proper par with the 
semiclassical method. As will be shown in this paper, 
this results in everything "falling into place." In partic- 
ular, the quantum potential — far from being singular in 
the vicinity of nodes — is well-behaved everywhere, and 
in fact, becomes vanishingly small in the large action 
limit, exactly in accord with the correspondence princi- 
ple. The same can be said for the quantum trajectories, 
which are no longer stationary, and approach the cor- 
responding semiclassical trajectories in the large action 
limit (within the classically allowed region of space) . This 
implies the somewhat counterintuitive result that quan- 
tum trajectories must be well-behaved when the number 
of nodes is large, for this signifies the large- action limit. 
In any event, the two-term Bohmian approach provides 
us with the "best of both worlds," i.e. the well-behaved 
trajectories of semiclassical mechanics, together with the 
singlevaluedness and lack of caustics and phase disconti- 
nuities that characterize Bohmian mechanics. More gen- 
erally than for just the stationary states in 1 DOF con- 
sidered here, it is anticipated that a multi-term Bohmian 
implementation will go a long way towards alleviating 
the node problem. 
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II. BACKGROUND 

A. Unipolar ansatz 

Let ^(x) be any normalized wavefunction in the single 
DOF, X. Being a complex function, '^/{x) can be de- 
composed into two real functions, Rb{x) and SBix)/h, 
representing the amplitude and phase, respectively, as 
follows: 



iSBix)/^. 



(1) 



Equation ^ is the celebrated Madelung-Bohm 
ansatz, ^^'^^'^^ which we term the "unipolar ansatz,"— 
as it consists of a single term only. The function Sb{x) 
has units and interpretation of action. If time evolution 
is considered, then Sb{x) plays the role of Hamilton's 
principle function in classical mechanics,— which there- 
fore properly depends on t as well as on x. However, for 
stationary states in a time-independent context, Sb{x) is 
analogous to Hamilton's characteristic function, W{x), 
which is time- independent. Both interpretations will be 
found to be important for the present approach. 

The above decomposition is essentially unique if Rb{x) 
is nonnegative throughout; however, it leads to dis- 
continuities in Sb{x), and cusps in Rb{x), if "^{x) 
has nodes. Despite these drawbacks, Eq. ^ is the 
decomposition generally utilized in standard Bohmian 
treatments (^iS^i^^ thus motivating use of the "i?" sub- 
script. It will be shown in this paper — evidently for the 
first time — that this convention in and of itself gives rise 
to certain node-related numerical difhculties that would 
otherwise not arise (Sec. Ill C|) . 

Accordingly, for the present work, we presume ampli- 
tude functions that change sign when passing through 
nodes. To avoid confusion with the standard Bohmian 
decomposition, we use unsubscripted quantities to repre- 
sent the present unipolar ansatz decomposition, i.e. 



(2) 



where Rb(x) — Sb{x) — S{x) mod (ttTi), and 

the new decomposition functions, R{x) and S{x), are 
smooth and continuous throughout the entire coordinate 
range; for the latter reason, S{x) is deemed a better 
analog for Hamilton's functions than is Sb{x). In any 
event, throughout this work, when discussing the unipo- 
lar ansatz, we are referring to Eq. ([2]), unless explicitly 
stated otherwise. Although the function S{x) is unique 
modulo 27r?i, physically, it is well-defined only up to the 
addition of an arbitrary constant, which introduces an 
arbitrary but immaterial phase factor into '^{x). 

Equation ([2]) is the starting point of both quantum tra- 
jectory and semiclassical methods. Presuming a quantum 
Hamiltonian of the form 



H 



V 

2m 



V{x), 



(3) 



the general time evolution equations for the Eq. ([2|) 
unipolar decomposition functions are 



R = -^(2R'S' 
2m ^ 

S = -V+^ 
2m 



RS") 



(4) 
(5) 



where R' — d [R{x)] /dx, R = d [R{x)] /dt, etc., and the 
coordinate dependences have been suppressed to save 
space. Equation ^ is the continuity equation, essen- 
tially stating that probability is conserved. Note that 
this equation is independent of the particular system po- 
tential, V{x). 

Equation ([5]) is the quantum analog of the Hamilton- 
Jacobi equation, 52 which does depend on the particular 
V{x). It is in the treatment of this equation that the 
semiclassical and Bohmian theories part company. The 
former ignores the first term in the square brackets, giv- 
ing rise to the standard classical Hamilton- Jacobi equa- 
tion. The latter regards the first term as the "quantum 
potential," 



Qix) 



2m 



Rr_ 

R 



(6) 



which is added to the true potential, to form the modified 
potential, U{x). Apart from the substitution V{x) 
U{x), the dynamical laws for the two approaches are 
identical. 

In particular, in both cases, the momentum P{x) is 
related to the action via 



P{x) = S'ix) 



(7) 



The set of points {x, P{x)} constitute a one-dimensional 
subspace of the 1 DOF phase space known as a "La- 
grangian manifold" (LM).-'-'^ J^i-'-^i-'-^ This terminology is 
generally used in a classical or semiclassical context only; 
however, we find it convenient to apply it in the exact 
quantum case as well. 



B. Stationary states, and the bipolar ansatz 

If ^'(x) is presumed to be a stationary state — i.e., an 
eigenstate of the Hamiltonian, _ff , with energy E — then 
R = 0, and S = ~E. The first result, together with 
Eq. (jH), is consistent with the well-known property that 
for standing waves, the phase is constant over x. Without 
loss of generality, we may take "^{x) to be real, so that 
S{x) — and R{x) — '^{x). The second result, together 
with Eq. yields the Quantum Stationary Hamilton- 
Jacobi Equation (QSHJE). By rearranging Eq. ([S]), and 
making use of Eqs. ^ and ([7]), we obtain: 



P\x) 



2m [E - 
2m \E 



V{x) - Q(x) 
U{x)] 



(8) 
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An important connection between the semiclassical and 
Bohmian theories is suggested by Eq. ([5]) — namely, the 
semiclassical approximation is accurate when the quan- 
tum potential is small. 

Let us now consider the semiclassical approximation 
proper — i.e., Eq. ((8]) with no Q{x) term. The resultant 
algebraic equation has two solutions for p{x), i.e. p{x) = 
±^y2m[E — V{x)] — where lower case 'p' is now used, for 
reasons that will be explained shortly. These two solu- 
tions correspond to the positive and negative momentum 
"sheets" of the LM in phase space. The two sheets are 
joined together at the classical turning points, Xmin and 
a^max, to form a single LM in phase space, along the classi- 
cal Hamiltonian contour, H{x,p) = p^/2m + V{x) — E. 
(Turning points are the caustics for stationary states). 
Thus, the function p{x) is double-valued over the classi- 
cally allowed region, Xmin < x < Xmax, and zero- valued 
everywhere else. 

It is illuminating to compare the semiclassical situa- 
tion, as described above, with the Bohm prescription. 
From the correspondence principle, one expects Q{x) to 
be small in the large action limit — i.e., the limit in which 
the excitation numbers, n, become large. This in turn 
would imply, via Eq. ([8]), that the semiclassical and ex- 
act quantum LMs would resemble each other in the large 
action limit. The actual situation is completely differ- 
ent, however. Firstly, Eqs. ([2]) and ([7]) imply that the 
quantum P{x) is single-valued everywhere — rather than 
zero- or double- valued, like the semiclassical p{x). Sec- 
ondly, for stationary states, S{x) = P{x) = 0, imply- 
ing that the quantum LM is the real axis of the phase 
space plane, which in no way resembles the semiclassical, 
Hamiltonian-contour LM. Thirdly, quantum trajectories 
are stationary over time, whereas semiclassical trajecto- 
ries evolve around their LMs (Sec. IIII D|) . 

The origin of these seemingly profound qualitative dif- 
ferences is deceptively simple: it is due to the fact that 
the semiclassical approximation does not actually incor- 
porate the unipolar ansatz of Eq. ([2]). Instead, a bipolar 
ansatz is used for the total wavefunction, consisting of 
two terms rather than one: 

■^{x) = e*V(a;)e*'*(^)/'' + e-'*r-(x)e~*"(^)/^ 

= +{x) + ^{x) (9) 

In Eq. ([9|) above, the constant 5 specifies the relative 
phase between the two components, '^±{x). It has been 
singled out from s{x) in order to clarify certain issues per- 
taining to square- integrability (Sec. IIIIB2|) . Although 
"^{x) itself is real, the '^±{x) are both complex, and con- 
jugate to each other. The stationary "standing wave" is 
therefore obtained in practice as the linear superposition 
of two "traveling waves," moving in opposite directions. 

The standard Bohmian prescription, being unipolar, 
completely misses out on this elegant and useful aspect 
of the semiclassical approach, which gives rise to a qual- 
itatively very different kind of amplitude/action decom- 
position. On the other hand, if the bipolar ansatz of 
Eq. ^ is incorporated into the Bohmian theory, rather 



than Eq. ([2]), then it is indeed possible to reconcile these 
two approaches, in a manner consistent with the corre- 
spondence principle, as will be shown in Sec. IIIII Note, 
however, that an important difference between semiclas- 
sical and bipolar quantum LMs is already evident in 
Eq. jl]); namely that the bipolar momentum function 
p{x) — -iis' [x) is double- valued throughout the entire co- 
ordinate range. Thus, the two exact quantum LM sheets 
never join, but extend into the classically forbidden re- 
gions all the way to a; = ±oo. 

Throughout this paper, we use lower case to denote the 
bipolar ansatz functions, so as to distinguish these from 
the unipolar ansatz functions, for which upper case is 
used. For convenience, all bipolar functions are hereafter 
defined to be single- valued everywhere, by referring to the 
positive-momentum LM sheet only — e.g., p{x) = s'{x). 



C. Node issues 

The increased flexibility of the bipolar ansatz is ex- 
tremely useful vis-a-vis the treatment of nodes, for it 
allows for the direct representation of nodes as interfer- 
ence fringes arising naturally between the two traveling 
waves. This possibility is exploited to great effect in semi- 
classical methods, which manage to contrive (approxi- 
mate) bipolar amplitude functions r{x) that are com- 
pletely nodeless — no matter how many nodes are present 
in "^(x) itself. Thus, apart from discontinuities near 
turning points (associated with Maslov indice&i^i^iii!^) , 
the semiclassical r{x) is smooth and positive, and the 
semiclassical s(x) is smooth and monotonically increas- 
ing. Moreover, these decomposition functions tend to be 
very slowly varying, in relation to ^>{x) itself, particularly 
when the latter has many nodes. 

The above properties would of course also be benefi- 
cial for exact QTMs — which from a practical standpoint, 
is a primary reason why the bipolar ansatz ought to be 
considered within a Bohmian context. One difficulty is 
that the exact quantum bipolar decomposition of Eq. ([9]) 
is not unique, in the sense that the QSHJE of Eq. ([S]) has 
a two-parameter family of solutions i^°'^^ In particular, 
one trivial solution is r{x) — "^{x); s{x) = = 6, which 
simply reproduces the unipolar result. Clearly, this is not 
the solution that we want, i.e., one that exhibits semiclas- 
sical correspondence in the large action limit; obtaining 
the latter will be the focus of Sec. IIIII 

For the unipolar ansatz, it has been stated that nodes 
always give rise to infinities in Q{x), owing to the singular 
denominator in Eq. ([6]). However, this is only strictly true 
if the standard Bohmian ansatz of Eq. ([T]) is used, for 
which the Rb{x) > convention is employed. If instead, 
one adopts the Eq. ([2]) convention, so that R{x) smoothly 
changes sign as a node is traversed, then Eq. ^ need not 
always be singular at a node. 

In particular, Q{x) is newer singular when "^(x) is a sta- 
tionary eigenstate of H, provided V{x) is well-behaved 
everywhere. This is because the time-independent 
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Schrodinger equation guarantees that the nodes of ^{x) 
also happen to be inflection points. Using Eq. ^ more- 
over, it can be shown that Q{x) = E — V{x). Thus, for 
the stationary states considered in this paper, even the 
unipolar ansatz is not singular, contrary to what previ- 
ously has been widely considered to be the case. Even if 
the standard Bohmian ansatz is used, Qb{x) at a node 
does not exhibit a singularity per se, but is rather ill- 
defined, owing to the cusp in Rb{x); away from the nodal 
point, Qb{x) = Q{x). 

In any event, we find it useful and convenient to distin- 
guish between two types of nodes, depending on whether 
Q{x) is formally well-behaved ("type one" nodes) or sin- 
gular ("type two" nodes). From a numerical perspec- 
tive, even type one nodes will cause difficulties for stan- 
dard quantum trajectory calculations performed using 
the unipolar Bohmian ansatz. This is because the slight- 
est numerical error in the evaluation of the Eq. ^ ra- 
tio will result in instability near the nodes — even though 
formally, Q{x) does not diverge. In contrast, due to the 
smoothness and lack of nodes of the r{x) functions that 
arise in the bipolar decomposition, numerical evaluation 
of the corresponding bipolar quantum potentials, q{x), 
causes no such instabilities for type one nodes. More 
general nodal implications of the bipolar ansatz will be 
discussed in greater detail in future publications. 



2. The classical probability distribution does not 
change over time. 

3. The classical flow for either LM sheet maintains 
invariant flux over all x, with the flux value for the 
two sheets being equal and opposite. 

4. The area enclosed within the LM, i.e. the enclosed 
action, J, is given by J = 2'Kh{n A- 1/2), where n is 
the number of nodes. 

5. For a normalized distribution, the absolute value of 
the invariant flux equals the inverse of the period 
of the trajectory, i.e. F = ujII-k. 

6. The median of the enclosed action, xq, satisfies 

7. All trajectories move along the LM. 

By "translating" these properties appropriately into the 
exact quantum context, we will be able to define an essen- 
tially unique bipolar decomposition of the Eq. ([5]) form. 

B. Basic properties: (1)— (4) 

1. invariant flux 



III. BIPOLAR DECOMPOSITION FOR 
STATIONARY STATES 

A. Semiclassical properties 

In this section, we derive a unique bipolar decom- 
position of the Eq. © form, for any given station- 
ary wavepacket, vE'(x), which satisfies semiclassical cor- 
respondence in the large action limit. In general, the 
Eq. ^ decomposition is nonunique. The semiclassical 
solution, however, is essentially unique fSec. Ill B|) . We 
will therefore use the latter as a guide, for selecting the 
particular quantum bipolar decomposition which most 
closely resembles the semiclassical solution. 

Note that certain assumptions have already entered 
into the form of Eq. ([5]), which is clearly more constrained 
than a completely general bipolar decomposition of ^(a;) 
into two arbitrary components. In particular, we have 
presumed ^{x) to be a superposition of equal and oppo- 
site traveling waves — a natural assumption, completely 
analogous to the semiclassical situation. Be that as it 
may, there is still an enormous number of ways in which 
Eq. ^ may be be realized for a given real wavepacket 
^(x), and so the decomposition is still far from being 
unique. 

To help narrow the field, we first summarize some of 
the additional properties of semiclassical eigenstates in 
1 DOE, which we will attempt to emulate in the exact 
quantum decomposition: 

1. The LM itself does not change over time. 



Properties (1) and (2) above are the most fundamental, 
and will be considered first. For a particular Eq. ^ de- 
composition, the corresponding positive-momentum LM 
sheet is given by p{x) = s'{x). Property (1) states that 
'p{x) = 0, which in turn implies s{x) = const. Property 
(2) implies r{x) = 0. Together, these properties imply 
that the ^±(a;) components of '^'{x) must be stationary 
eigenstates of H in their own right. Equation ^ then 
implies that the eigenvalues for the two components must 
both be equal to E. 

Since the '^±{x) are stationary, the quantum mechan- 
ical flux associated with each of these components, i.e. 



j±{x) — 

2im 



= ± 



p{x) 



d'^±{x) d-iilix) 



dx 



dx 



^±{x) 



(10) 



is independent of x, with equal and opposite constant 
values, ±F. Note that Eq. (fTO)) demonstrates that the 
quantum analog of property (3) is also satisfied — a nec- 
essary consequence of properties (1) and (2). We call this 
the "invariant flux" property. 

For F = 0, ^'+(x) = 4'_(x) (X ^'(x), which reproduces 
the unipolar ansatz. We shall therefore hereafter restrict 
consideration to the F > Q case, for which the invariant 
flux property, and Eq. (jlO[) . provide a speciflcation for r, 
in terms of s' and the constant, F: 



-(x) 



' mF 
s'{x) 



(11) 
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Note that if r{x) > for all x, then Eq. ([TT|l implies that 
p{x) > for all X — a desirable property for the positive 
momentum solution, also satisfied by the semiclassical 
solution. This would also imply that s{x) is monotoni- 
cally increasing. Accordingly, the r{x) > condition is 
adopted. 

For F > 0, the two 5'±(x) components are linearly 
independent. This implies that at least one of the two 
must be non-L^. In fact, being complex conjugates of 
each other, both solutions must be non-L^. Moreover, it 
can be shown that 4'±(x) diverges as a; ±oo.— This 
is due to the fact that p{x) = s'{x) ^ as a; — > ±oo in 
order that the enclosed action be finite fSec. lIIIB2|) : but 
this implies via Eq. (fTT|) that r{x) diverges. 

It is instructive to rederive the invariant flux property 
in another manner. By the superposition principle, the 
time evolution of ^>{x) can be obtained by propagating 
each of the ^'±(a;) components separately in time. Since 
the two components are stationary eigenstates of H in 
their own right, the time-evolving \E'-(-(x)'s must each in- 
dependently satisfy r = and s = —E. The former, 
applied to a lower case version of Eq. (jl]), is equivalent 
to the spatial derivative of Eq. (fTO|) . 

2. quantization 

We now address the exact quantum analog of property 
(4), the quantization condition. In a proper semiclassi- 
cal treatment, this half-integer condition on the enclosed 
actiouf^i^iS J = 2TTh{n + 1/2), must be supplemented 
by the discontinuous jumps in phase that occur as one 
traverses a turning point, from one LM sheet to another. 
In a certain sense, these jumps account for the fact that 
the WKB solutions do not incorporate the portion of 
the true wavefunction that tunnels into the forbidden 
region — which contribute an additional one half quanta 
of action, over the course of one complete circuit around 
the LM,-'^^ When this discontinuous contribution is prop- 
erly added to the usual enclosed action contribution, one 
obtains an integer quantization condition for the total ac- 
tion, Jtot = 27r?i(n-|-l), even within a purely semiclassical 
context. 

In the quantum case, there is no distinction between 
classically allowed and forbidden regions; one travels 
smoothly from one to the other, over the entire position 
space. Since the two LM sheets are symmetrically placed 
in phase space about the real axis (Fig. [T]), the area en- 
closed between them is clearly twice the change in the 
action function. As = J^""^'' p{x)dx, as one travels from 
a^min = ^oo to Xniax = +00. From the above description, 
we expect this change in action to be nhin -1-1), where 
n is the number of nodes. An integer quantization con- 
dition is therefore expected to hold for the exact bipolar 
quantum decomposition. This is indeed correct, as has 
been shown previously .^'^'^'^i^^ 

A sufficient, though certainly not necessary (see be- 
low), condition for achieving integer quantization of the 



(a) 



T 1 < 1 T 




position X 

FIG. 1: Bipolar Lagrangian manifolds (LMs) for the harmonic 
oscillator ground state, vl'(a;) = exp(— a;^/2)/7r^''^, for three 
different flux values, F {xq = in each case): (a) F = l/2n « 
0.159, the semiclassical value; (b) F = 1/A-/t^ ^ 0.141; (c) 
F — l/2\/27r ~ 0.199. Solid curves indicate quantum LMs; 
dotted curves indicate semiclassical LMs. The former en- 
close an area J = 2nh; the latter J — nh. In (a), the small 
open/filled circles represent semiclassical/quantum trajectory 
locations at times t = kTr/2, for integer k. 

quantum action is that ^/{x) be L^, which requires that 
^(±cx)) = 0. For convenience, we adopt the convention 
that s{-oo) = -nh{n + l)/2 in Eq. The condition 
^'(— oo) = then determines a value for S — i.e., 6 = for 
n even, and S = tt/2 for n odd. This yields the following: 

^(x) = / cos[s(x)/n] for n even; 
^ ' 1 2r(a;) sin[s(x)/?i] for n odd. 

The somewhat awkward distinction between the even and 
odd n cases is due to our choice of boundary condition 
for s(x)\ the reasons for this seemingly perverse choice 
will be made clear shortly. 

Since ^'(-l-oo) must also be zero if *(x) is L^, Eq. IT^ 
implies that s(-l-oo) = ■nfiiri! -\- 1)/2, where n' is an inte- 
ger. Since rix) > everywhere, all nodes of '^{x), for n 
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even/odd, must occur at x values for which s{x)/fi is an 
odd/even multiple of 7r/2. This, in turn, requires n' = n, 
from which is obtained 

J = 2As = 2s(+oo) - 2s(-oo) = 2T:h{n + 1), (13) 

i.e., the integer quantization condition. Note that s{x) is 
independent of the normalization of '^{x). However, the 
amplitude function r{x) is not; thus, if '^'{x) is presumed 
normalized to unity, then Eq. p2|) defines the normaliza- 
tion convention for r{x), and for '^±{x). 

The specific value of the constant, S — as given above 
for even/odd n — is not arbitrary, but is required in order 
to ensure that the resultant ^'(x) be L^. If a different 
value were chosen, then a different, non-L^ solution to 
the Schrodinger equation would be obtained. This situ- 
ation is in stark contrast to the unipolar case (for which 
the constant simply changes the overall phase) , and is due 
to the fact that the bipolar 6 represents a relative phase 
shift. In any event, we can regard 5 as a parameter that 
is used to specify a particular solution of the Schrodinger 
equation at the energy E. In fact, all such real-valued 
solutions (apart from an immaterial scaling factor) may 
be obtained by varying the value of 5 in Eq. ©. This 
is due to the fact that the '^±{x) are linearly indepen- 
dent. In any event, an important consequence is that all 
Schrodinger solutions are described by the exact same 
action function, s(x), which is unaffected by the value 
of 5. Among other things, this implies that the integer 
quantization rule applies to all non-L^ solutions, as well 
as to the solution. 



C. QSHJE properties: (5) and (6) 

1. introduction 

In this section, we continue the approach introduced at 
the end of Sec. lIII BTI where the time evolution equations 
are applied to the two '^±{x) components separately. 
The r = equation was seen to yield the same invariant 
flux relation, for each 5'±(x) component — i.e., Eq. pT|) . 
The s = —E equation, as applied separately to the two 
components, also yields identical results. One of the im- 
portant ramifications of this is that the Eq. ^ ansatz is 
preserved over time, at least for stationary eigenstates. 

In any event, since '^±{x) are solutions to the 
Schrodinger equation, the s = E equation must result 
in a lower-case version of the QSHJE [Eq. ([5])], for which 
the bipolar quantum potential, q{x)^ is defined via lower- 
case Eq. (ini). We shall rewrite this QSHJE by expressing 
q{x) directly in terms of p{x), using lower-case Eqs. ^ 
and dZl), and Eq. pl|)i^ 



= 2vi{E -V)-h^ 



1 p1 

2 \p 



I p1 

4 \p 



solution, Eq. p4|) is a second-order differential equation 
in p{x) , with a two-parameter family of different solutions 
to choose from. 

Note that Eq. ([H]) applies to all Schrodinger equation 
solutions, i.e. the and the non-L^ solutions, both. 
Since we are interested only in the former, and since 
there is a one-parameter family of Schrodinger solutions 
in all, one might expect that the specification of the 
solution would determine the value of one of the two pa- 
rameters. This is not correct however, as demonstrated 
earlier (Sec. lIIIB 2p . Thus, even for the solution alone, 
bipolar decomposition gives rise to two variable parame- 
ters via Eq. (fll)) . 



2. defining the two parameters 

What are the two parameters, and how should their 
values be chosen? To determine what the two parame- 
ters are, it is convenient to combine Eqs. pT|) and (fT^ 
together, to obtain a formula for ^i{x) directly in terms 
of s{x). The result is: 



^(2;) = I \ /^TTT'F/s'{x) cos[s{x)/h] for n even; ^^^^ 
[ yjAmF/s'{x) sin[s(2;)/?i] for n odd. 

Equation (|15|) is a first-order differential equation for 
s{x); the general solution is easily found to be 



, 4mF r dx' 

tan(s/?i) ~ — - — / 



- coi{s/h) 



5-2(2;') 
AmF r dx' 



for n even; 



for n odd. (16) 



(14) 



For nodeless wavepackets {n — Q), 5-(x) > everywhere, 
and the integrand of Eq. (fT6|) has no singularities. When 
n > 0, Eq. (|16p is still correct, but requires careful branch 
selection, to ensure that the final s{x) curve is continu- 
ous throughout the coordinate range. Note that 5- (a;) is 
presumed to be the Schrodinger solution. 

Equation (|16p provides an explicit recipe for obtain- 
ing the Eq. (|12p decomposition. The two parameters can 
thus be taken as: (1) the flux parameter, F; (2) the inte- 
gration limit parameter, xq. Note that s(xo) = 0; conse- 
quently, xq may also be interpreted as the median of the 
action, as per Sec. IIII Al By varying the two parameters 
F and xq in Eq. (|16p . different bipolar decompositions 
may be achieved. These correspond to different affine 
transformations of each other, in the sense that varying 
F is equivalent to rescaling the right-hand-side (RHS) of 
Eq. (jl6p . whereas varying a;o is equivalent to adding a 
constant to the RHS. 



choosing parameter values 



Whereas the semiclassical approximation [obtained by ig- 
noring the q{x) term in brackets] has a unique p{x) > 



For a given 5'(x), the various bipolar LMs that can be 
constructed via Eq. pop vary significantly with respect 
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to F and xq fSec. lIV|) . and so a general procedure for ob- 
taining reasonable parameter values must be provided. 
At present, the best approach seems to be to touch base 
once again with the semiclassical properties — in partic- 
ular, property (5) for determining the appropriate value 
of F, and property (6) for determining the appropriate 
value of Xq. 

Semiclassically, the flux for a normalized distribution 
is given by = a;/27r, where lo is the classical angular 
frequency for the appropriate semiclassical trajectory — 
i.e., the (uniform) rate at which the angle variable of the 
action/angle pair changes. The corresponding quantum 
trajectory is not that of a bound state, and so it is not 
possible to assign an oj value to it (Sec. IIII Pp . On the 
other hand, the Eq. ()12|1 normalization convention al- 
lows us to determine a unique flux value for the quantum 
trajectory, which is all that is required. By setting the 
quantum flux value equal to the semiclassical value, it is 
anticipated that the resultant quantum LMs will closely 
resemble the semiclassical LMs, as desired. 

As for Xq, the median of the enclosed action: this can 
be can be regarded as the exact middle of the wavepacket 
in a certain sense; semiclassically, Xq is the classically 
allowed configuration that is furthest from both of the 
turning points, vis-a-vis the action. Consequently, we 
expect the greatest agreement of semiclassical and quan- 
tum LMs — i.e. the smallest q{x) values — in the vicinity 
of the semiclassical xq. This can be achieved by allow- 
ing the quantum a;o to coincide with the semiclassical 
value — i.e., the latter is chosen to be the location where 
s = 0. 

In the quantum bipolar decomposition scheme — even 
for fixed F — one is otherwise free to place the action me- 
dian, Xq, essentially anywhere along the position axis. 
The ramifications are particularly illuminating when 
V{x) is even. For such potentials, only the choice a;o = 
gives rise to quantum bipolar decomposition functions 
that are even or odd in x, thereby respecting the physical 
symmetry of the system, and of '^{x) itself. This choice 
for Xq is also consistent with the median action criterion. 
Presumably, it would be unphysical to consider any of the 
asymmetrical decompositions; nevertheless, it is interest- 
ing to note that one can generate asymmetrical bipolar 
decompositions that give rise to the symmetrical 4'(x), 
simply by shifting xq away from the origin. This has 
been verified via explicit construction for the harmonic 
oscillator ground state. 



D. Quantum trajectories: property (7) 

We now address the issue of the quantum trajecto- 
ries themselves, related to property (7). Semiclassically, 
over the course of time, the bound state trajectories 
simply move around and around the Hamiltonian con- 
tour LMs, which do not themselves change [property 
(1)]. In a conventional unipolar quantum treatment, 
the initial LM — specifying the initial conditions for the 



ensemble of trajectories — is just the real axis, i.e. the 
"curve" P{x) = 0, since '^/{x) is real. The quantum 
trajectories evolve under the unipolar modified poten- 
tial, U{x) = V{x) + Q{x), which by Eq. must be 
the constant function U{x) ^ E [even if there arc nodes 
fSec. IIICp ]. Consequently, P = -U'{x) = 0, and so the 
unipolar quantum trajectories do not move at all over 
time. 

In contrast, the bipolar quantum trajectories are not 
stationary, but move along the positive and negative mo- 
mentum LM sheets. This is true because p{x) > 0, 
and because the bipolar LMs themselves do not change 
over time, thus verifying property (7). Moreover, pro- 
vided the bipolar quantum potential q{x) is small in 
the classically allowed region, then the bipolar quan- 
tum trajectories must resemble the semiclassical trajec- 
tories within this region, since the LMs are similar, and 
p = —u'{x) w —V'{x). Of course, the bipolar quantum 
trajectories do not change their direction at the classical 
turning points, moving between positive and negative- 
momentum LM sheets, like classical trajectories. Instead, 
all quantum trajectories on say, the positive momentum 
LM sheet, continue to head to the right for all time. Once 
these trajectories enter the classically forbidden region, 
however, their speed decreases very suddenly. 

It is worth discussing the very different role played by 
the quantum potential in the unipolar ansatz, versus that 
of the bipolar ansatz with the specific decomposition sug- 
gested here (i.e., parameter choices of Sec. IIII C3^ . In the 
unipolar case, Q{x) serves to counteract the true poten- 
tial everywhere; thus, Q{x) is not generally small. In 
contrast, the bipolar quantum potential, q{x), can be 
regarded as the 0{h^) correction to the semiclassical 
approximation — in the truest correspondence-principle 
sense of lower-case Eq. (|8]). The value of q{x) is therefore 
small in the appropriate semiclassical limits — i.e., in the 
classically allowed region far from turning points, and in 
the limit of large action, when n becomes large. Near 
the turning points, \qix)\ increases substantially, so as to 
ensure that all trajectories keep moving past the classical 
turning point without changing direction. This increase 
continues well into the classically forbidden region, where 
curiously, q{x) approaches Q{x) — i.e. it effectively can- 
cels out the true potential. Consequently, the bipolar 
modified potential, u{x) = V{x) + q{x), resembles the 
true potential in the classically allowed region, and the 
unipolar modified potential, U{x) — E, in the asymptotic 
regions. 

The above discussion hinges on the assumption that 
the semiclassical and bipolar quantum LMs become ar- 
bitrarily close in the appropriate semiclassical limits de- 
scribed above. We can justify this expectation as fol- 
lows. First, the semiclassical approximation is known to 
become arbitrarily accurate in these limits; each of the 
two semiclassical traveling wave components must there- 
fore approach some particular corresponding pair of ex- 
act quantum solutions, '^±{x), arbitrarily closely. The 
latter must therefore have the same characteristics, vis- 



9 



a-vis action, trajectories, and flux, as do the semiclassi- 
cal approximations, in the appropriate limits. Therefore, 
by choosing the available parameters for the quantum 
solutions (i.e. F and xq) so as to match the semiclas- 
sical approximations, the correspondence principle must 
be satisfied. 



E. Stationary non-eigenstates 

Although the primary interest of this paper is bound, 
stationary eigenstates of the Hamiltonian H of Eq. ([3]), 
our ultimate interest is wavepackets that evolve dynami- 
cally over time. As a first step in this direction, we gener- 
alize the previous discussion to include wavepackets that 
are only momentarily "stationary." In other words, the 
initial wavepacket ^(x) is real, but otherwise arbitrary, 
i.e., not presumed to be an eigenstate of H. This re- 
sults in r = 0, but only instantaneously, at the initial 
time. We shall call such a wavepacket a "stationary non- 
eigenstate." 

To what extent can the bipolar decomposition scheme 
be applied to stationary non-eigenstate wavepackets? 
The question is relevant, because it is only necessary 
to specify the bipolar decomposition at a single point 
in time, in order to propagate the two vE'±(x) compo- 
nents independently, over all time. Our approach shall 
be to regard '^{x) as the eigenstate of some Hermitian, 
Hamiltonian-like operator, Hq, which without loss of gen- 
erality, may be taken to be of the form 



2m 



Voix) 



(17) 



If \I'(a;) is known, it is a trivial matter to obtain Vo(a:) by 
solving the Schrodinger equation in reverse, i.e. 



[Vo{x)-Eo] 



2m 



4(" 



(18) 



The Hq so obtained can then be used to generate the 
appropriate F and xq values semiclassically. 

As a completely general procedure, this approach has 
one unavoidable flaw. If \l/(x) has type two nodes, then 
Vo{x) will have singularities at the nodes, which is un- 
desirable. In such cases, since the precise values of the 
F and xq parameters may not matter all that much in 
numerical practice, one should simply choose "reason- 
able" values by comparison with known cases — e.g., for 
n — 0, one could use the parameters of a Gaussian with 
the same center and standard deviation as \E'(a:). On 
the other hand, almost all initial wavepackets used in 
chemical physics applications correspond [via Eq. (jlSp ] 
to potentials Vo{x) that are well-behaved. 



IV. RESULTS: HARMONIC OSCILLATOR 
EIGENSTATES 

As a classic benchmark example, we now work out an- 
alytic solutions for the harmonic oscillator (HO) eigen- 
states, i.e. V{x) = kx^ 12. This example is particularly 
important, as the ground state provides the proper de- 
composition for Gaussian wavepackets, which are used 
very frequently in time-dependent studies. We shall also 
find the excited harmonic oscillator states to be quite use- 
ful, particular with respect to investigations regarding 
nodes and interference. The normalized n'th harmonic 
oscillator eigenstate is given by 

1/4 



Fir, 



'/2h 



(19) 



where Hn{) is the n'th Hermite polynomial, and uj 
y^k/m. 



Ground state 



We start with the ground state, ^'0(2^) = 
(mw/n7r)i/4e-™"^'/2'*. Application of Eq. ^ yields 



s{x) — ?iarctan-| F ( — ) x 



rfi 



— crfi 



-xo 



(20) 



as the generic, F- and xo-dependent solution. This gives 
rise via Eq. pTjl to 



r{x) 



'1\ rmuj\^/'^ 
2) 



[ ^ ' ' ^~mujx''/2h ^ 



1 + F^ 



erfi 



2n 



(21) 



erfi 



-Xo 



2 \ 1/2 



As per Sec, nil C3[ the appropriate value of F is clearly 
F = w/27r. The appropriate value of xq, whether from 
symmetry considerations, or the more general median ac- 
tion criterion, is clearly xq = 0. With these choices for 
the parameter values, we obtain the simpler result 



s{x) 
r(x) 



h arctan 



erfi 



1\ /■muj\ 



1/4 



2 J \h-K J 

1 + CT& 



C X 



1/2 



(22) 



(23) 
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-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 
position X 

FIG. 2: Amplitude functions for tiie harmonic oscillator 
ground state, "^{x) = exp(— a;^/2)/7r^''*. Dashed curve: 
unipolar amplitude, R(x) = ^{x). Solid curve: bipolar am- 
plitude, r{x), for semiclassical parameter values, F = l/27r 
and So = 0. 



For convenience, we choose units such that h ^ m — k = 
uj = 1. In these units, F = l/27r, and Xq = 0. In these 
units, and for these parameter values, all of the relevant 
bipolar decomposition functions are as follows: 

s{x) — arctan [erfi (x)] 

r{x) = Q)(iy^%-V^[l + erfi^(.)]^/^ 



' ^ ' 0F[l + erfi2(x)] 

q{x) = - - — ^ (24) 

2 2 7r[l-Herfi2(x)]' 

All of the functions in Eq. ((24|) are smooth, slowly vary- 
ing, and monotonic in |a;|. The LM is an elegant "eye- 
shaped" curve [specified by the p{x) equation above] that 
deviates smoothly, and positively, from the circular semi- 
classical LM, with the point of closest approach being 
xo = 0. All of these features are as predicted in Sec. IIIIl 
and would not have been satisfied if substantially differ- 
ent parameter values were used. A plot of the semiclas- 
sical and bipolar quantum LMs is presented in Fig. [U 
for F — l/27r and other values (but all with xq — 0). 
Whereas some of these other plots have the qualitatively 
correct behavior, it is very clear that the F — \/2ti curve 
is the smoothest, most "correct" choice — especially vis-a- 
vis comparison with the corresponding semiclassical LM. 

Figure [2] is a comparison between the unipolar and 
bipolar amplitude fimctions — i.e., R{x) and r{x), respec- 
tively. As is clear from the figure, these two types of 
amplitude behave completely differently. In particular, 
whereas R{x) decreases quickly as x ^ ±oo, r{x) in- 
creases as one moves away from the origin, and actu- 
ally diverges in the x — > ±oo limits, as predicted in 
Sec, nil BT] Clearly, the 5'±(x) are non-L^ solutions. 



(a) 



T |—p 1 1 1 1 1 r-p r 

0.10 T-~, 1 1 h T—, 




position X 

FIG. 3: Bipolar modified potentials, u{x), and true poten- 
tial, V{x), for two different harmonic oscillator eigenstates: 
(a) ground state, n = 0; (b) tenth excited state, n = 10. 
Solid curves indicate u{x); dotted curves indicate V{x). In 
(a), u{x) approaches i?„=o = 1/2 as a: — > ±oo, but resembles 
V{x) in the classically allowed region, \x\ < 1 [inset shows 
u{x) — «(0) vs. V(a;)]. In (b), u{x) and V{x) are virtually in- 
distinguishable over the (classically allowed) coordinate range 
indicated. 



In Fig. [3fa), we present a comparison of the actual 
and bipolar modified potentials — i.e., V{x), and u{x) — 
V{x) + q{x). The two potentials resemble each other 
in the classically allowed region, away from the turning 
points at X = ±1. As one approaches the turning points, 
the difference q{x) increases markedly. In the classically 
forbidden region, u(x) ceases to emulate the true poten- 
tial, and in the asymptotic limits, approaches the unipo- 
lar U{x) constant value oi E = 1/2. All of this is as 
predicted in Sec. IIIIDI 

We also performed trajectory calculations. In particu- 
lar, a single trajectory was propagated over a very long 
period of time, using what Wyatt has called the "analyt- 
ical approach." In this scheme, the modified force [i.e. 
— u'(x)] is computed analytically, but the trajectory itself 
is propagated numerically. We found first of all that this 
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numerical propagation scheme was extremely stable, as 
demonstrated by the fact that the numerical trajectory 
did not deviate appreciably from the LM at any point 
in time. Phase space values for the trajectory at vari- 
ous times are indicated as small circles in Fig. [ija) , from 
which it is also clear that quantum trajectories corre- 
spond fairly well to the semiclassical trajectories in the 
classically allowed region. 

In the classically forbidden regions, trajectories slow 
down very quickly, as predicted. This is evidenced by the 
pile-up of trajectory points that ensues in these regions 
[Fig. [H^a)]. Formally, however, the trajectories do not 
actually reach zero momentum until x —> oo. They are 
thus analogous to classical trajectories for a system that 
has just enough energy for dissociation. This fact is also 
reflected in the asymptotic behavior of u{x) as discussed 
above. 

We now briefly address the issue of trajectory "pile-up" 
in the classically forbidden regions, which is an important 
concern for numerical calculations. Although the bipolar 
ansatz has the effect of placing more trajectories in re- 
gions of space where the actual probability is small, this 
situation is numerically agreeable for two reasons: (1) 
more accuracy is needed in these regions, because "^{x) 
itself is effectively obtained as the difference between two 
large numbers; (2) a simple "recycling" scheme can be in- 
troduced to reduce the number of trajectories to a min- 
imum. These issues will be discussed in great detail in 
future publications. 



B. Excited states 

The correspondence between the semiclassical and 
quantum LMs for the harmonic oscillator ground state is 
only fairly good, but one ought to recall that the action 
is minimal in this case. A real test of the correspondence 
principle requires a detailed investigation of the LM be- 
havior in the large action limit. This in turn, requires 
that the bipolar decomposition be performed for the ex- 
cited harmonic oscillator states. Using Eq. with 
h — m — k = uj = l, F — l/27r, and xq = 0, we have 
obtained analytical solutions for all n up to rt = 10. 

The general form of the bipolar action for the nth 
eigenstate [denoted Sn{x)] is as follows: 



Sn{x) 



arctan 



-arccot 



+ erfi (x) 
erfi (x) 



for n even; 
for n odd. 



(25) 

In Eq. (f25|) above, fn{x) is an (n — l)th-order odd/even 
polynomial, and gn{x) is an nth-order even/odd poly- 
nomial, for n even/odd. Explicit coefficient values for 
fn{x) and gn{x) are listed in Tables HI and HIl respectively 
(coefficients for larger n can be provided on request). 

The bipolar solutions for the excited states behave ex- 
actly as predicted. In particular r„(a;) > everywhere. 




-2.0 0.0 2.0 
position X 

FIG. 4: Bipolar Lagrangian manifolds (LMs) for three dif- 
ferent harmonic oscillator eigenstates; moving concentrically 
outward from the origin, these are n = 0, n = A, and n — 10, 
respectively. Solid curves indicate quantum LMs; dotted 
curves indicate semiclassical LMs. The correspondence prin- 
ciple is clearly satisfied with increasing n. 



and both rn{x) and Sn{x) are as smooth, slowly varying, 
and monotonic as for the ground state. In fact, all of the 
bipolar functions qualitatively resemble those for n — Q, 
except on a larger scale — as is also true of the semiclassi- 
cal functions. This is in sharp contrast to the behavior of 
the wavefunctions, ^'„(x), themselves, which gain nodes 
and rapid oscillations as n is increased. One very en- 
couraging aspect of the excited state bipolar functions is 
that nodal features are not evident anywhere. However, 
this requires choosing the correct branch of the Eq. p6|) 
solution at every point in position space, such that the 
resultant Sn{x) curve is not discontinuous across a node. 
This issue is discussed in more detail in Sec. IIV CI 

In any case, the basic goal of the bipolar ansatz 
has been achieved — i.e., to obtain a decomposition 
which, like the corresponding semiclassical approxima- 
tion, treats all nodes as interference between a superpo- 
sition of left and right traveling waves, '^±{x), which are 
themselves nodeless. Moreover, beyond achieving just 
this basic goal, we find that the correspondence princi- 
ple is satisfied in the large n limit. This is exemplified 
in Fig. m wherein the semiclassical and bipolar quantum 
LMs are presented for several harmonic oscillator states 
over the n range considered. 

From the figure, the quantum LMs are seen to enclose 
one-half quanta of area more than the semiclassical LMs, 
which manifests primarily in the forbidden regions near 
the turning points, as expected. In a relative sense, this 
discrepancy becomes decreasingly significant in the large 
n limit. Note that the quantum LMs completely enclose 
the semiclassical LMs, which was not anticipated earlier, 
but is certainly a desirable property. From lower-case 
Eq. ([8]) , this will only be satisfied if the bipolar quantum 
potential is negative everywhere. This has been observed 
for all examples considered thus far, using the appropri- 
ate semiclassical values for F and xq; however, most other 
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TABLE I: Coefficients, aj, for polynomials, fn{x) = X/"-o ^^^^^ all < n < 10. These polynomials appear in the numerator 

of the analytic formula for the bipolar action, s„(a;) [Eq. (I25|l ]. 
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TABLE II: Coefficients, bj, for polynomials, gn{x) = ^"^pfoja;-', for all < ti < 10. These polynomials appear in the 
denominator of the analytic formula for the bipolar action, Sn{x) [Eq. (|25|l ]. 
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parameter choices would not satisfy this property. 

Even more so than in the ground state case, the bipo- 
lar quantum potential for (x) is found to be very small 
throughout most of the classically allowed region. The 
magnitude of qn{x) decreases with increasing n, so that 
whereas go(0) = —.137 for the ground state, by n = 10 
we have gio(O) = —.012. Of course, the extent of the 
classical region is also larger with increasing n; thus by 
n = 10, we find wio(x) to be practically indistinguishable 
from the V{x) over the range |a;| < 3. The situation, 
depicted in Fig. [21Ib), can be regarded as another mani- 
festation of the correspondence principle. 

The correspondence principle also has important ram- 
ifications for trajectory calculations. In particular, not 
only are the bipolar quantum trajectories for large n 
smooth and well-behaved throughout, but in the classi- 
cally allowed region, they are virtually indistinguishable 
from classical trajectories. This has once again been ver- 
ified by performing analytical trajectory calculations for 
n = 10, which were found to be just as numerically sta- 
ble as for the ground state — despite the fact that 5*10 (x) 
itself has ten nodes. This bodes very well for obviating 
the node problem in general. 



C. Numerical issues 

Although the bipolar decomposition functions — once 
obtained — exhibit no special behavior in the vicinity of 
nodes, it turns out that nodes complicate the determi- 
nation of these functions somewhat, vis-a-vis implemen- 
tation of Eq. . To begin with, let us imagine that — 
as in the current harmonic oscillator case — an analyti- 
cal expression for the Eq. (flBl) integral is available. For 
the moment, we also take n to be even. Note that the 
left and right sides of Eq. must be infinite at the 
nodes. Thus, whereas the exact analytical expression can 
be used across the entire coordinate range, a new branch 
is encountered each time a node is traversed. The spe- 
cific branch of interest is specified by the condition of 
continuity for s„(a;), and by s„(a;o) = 0. 

A superficial difficulty is encountered for the odd n 
states, for which there is necessarily a node at xq. Strictly 
speaking, this implies that the Eq. (jl6p integration must 
be singular. To circumvent this difficulty, we express 
the RHS of Eq. (fTB]) as an indefinite analytical integral, 
plus an arbitrary constant, B. Note that since xq must 
lie at a node for odd n, xq can not serve as the sec- 
ond parameter, for singling out the particular solution of 
interest for Eq. We can, however, use B for this 
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position X 

FIG. 5: Bipolar Lagrangian manifolds (LMs) for the fourth 
excited eigenstate of the Morse oscillator system with twenty 
bound states total. Solid curves indicate quantum LMs; dot- 
ted curves indicate semiclassical LMs. Semiclassical values 
for F and a;o were used to specify the quantum solution, as 
per Sec. HVCl 



purpose. In particular, if Vix) is even, then only one 
value of B gives rise to the requisite odd s„(x) solution. 
More generally, i.e. for arbitrary V{x)^ we can still apply 
the oddness criterion locally. In other words, it is easy 
to show that B should in general be chosen such that 
Sn(a;o ~ e) = — Sn(a;o + e) in the limit of small e. This 
technique bears a resemblance to the Cauchy principal 
value method 

The above discussion presumes that an analytical in- 
tegral is available for the RHS of Eq. ITBl) . Generally 
speaking, this will not be the case, and wc must consider 
how to apply the above procedures when the integrations 
are performed numerically. Fortunately, this is straight- 
forward. The general procedure is to pick an arbitrary 
integration limit, Xi, lying in between each adjacent pair 
of nodes (where x = ±(X) are treated as "nodes" in this 
context). Equation is then applied to each interval 
separately, generating a smooth, numerically integrated 
function over the entire i'th interval, that is correct to 
within an additive constant, Bi. The individual Bi val- 
ues are then obtained, using the constraint s„(xo) — 0, 
and applying the Cauchy-like condition described above 
across each of the nodes separately. 

The numerical procedure described above has been 
applied succesfuUy to the Morse oscillator system, for 
which: (1) V^x) is not symmetrical, and; (2) the Eq. ifTB]) 
integrations must be performed numerically. The results 
will be presented in a future paper. We mention the 
Morse investigation in this paper simply to emphasize the 
fact the present method is in fact applicable in a much 
broader context than the analytical harmonic oscillator 
system considered here. Moreover, all of the conclusions 
drawn for the harmonic oscillator, regarding the corre- 
spondence principle and the like, evidently apply to more 
general systems. This is demonstrated in Fig. [Sj which 
depicts the semiclassical and bipolar quantum LMs for 
the n = 4 state of the Morse oscillator. 



V. COMPARISON WITH RELATED METHODS 

The present work does not constitute the first appli- 
cation of the bipolar ansatz in a Bohmian-like dynam- 
ical context. For two decades or so, a very interest- 
ing bipolar approach has been developed and advocated 
by E. R. Floyd-^i^i^i^i^i^ More recently, essentially 
the same technique was derived by Faraggi and Matone 
(FM)^i^i^ within a much broader context, and using 
a very different physical picture. It is somewhat remark- 
able that these two approaches give rise to exactly the 
same dynamical equations (an illuminating comparison is 
presented in Ref . [H and in Ref . [sil ) . Perhaps even more 
remarkable, however, is that the dynamical law used in 
Floydian/FM trajectory propagation is not equivalent to 
that of Bohmian mechanics. 

In this section, we compare and contrast the methods 
of Floyd, FM, and the present work. The various ap- 
proaches were developed independently, and so a brief 
discussion of the different philosophies is presented, as 
well as the mathematical similarities and differences. The 
comparison is particularly apt for the present paper, in 
that both the Floyd and FM theories are restricted to 
stationary states only. 

The starting point of the Floydian approach is the 
QSHJE in 1 DOF (multidimensional systems may be con- 
sidered, but only if there is separation of variables^-). 
As per Sec. [Ill this is natural enough, for a Bohmian- 
like theory applied to stationary states; however, there 
are some subtleties regarding the manner in which en- 
ergy is treated, that give rise to a non-Bohmian dy- 
namical law. The intriguing approach of FM begins 
not with the QSHJE, but with the fundamental pos- 
tulate that all systems are equivalent under coordinate 
transformations . ^^'^^ This is termed the "equivalence 
principle," in obvious analogy with general relativity, 
with which the FM approach shares many parallels. In- 
deed, one of the goals of FM is to reconcile quantum 
mechanics and general relativity, to which end it is natu- 
ral to focus on the action, which plays a key role in both 
physical theories. 

In any event, the FM approach relies on the existence 
of a particular coordinate transformation that reduces 
the system to that of a free particle i^i^ This approach 
implicitly relies on the fact that the Wigner-Weyl cor- 
respondence is not preserved under canonical/unitary 
transformations^^ — a feature that has been exploited to 
great effect in improvements to semiclassical theories, 
such as the Langcr modification , ^^i^''' although not al- 
ways with full comprehension that this is what was tak- 
ing place i^i^ 

Starting from the equivalence principle, FM derive the 
QSHJE, and demonstrate that energy quantization arises 
naturally from the condition that this postulate be sat- 
isfied. They find, moreover, that a unipolar ansatz is 
insufficient to achieve this, but rather, a bipolar ansatz 
of the Eq. ^ form is required. In contrast, Floyd begins 
with the QSHJE, and then presumes a bipolar ansatz 
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for which each of the two terms is a sohition. Floyd's 
motivations are evidently similar to the author's, in that 
he introduces the bipolar ansatz in order to !^°'^^ (1) ob- 
tain a Bohm-like formalism that is well-behaved at nodes; 
(2) obtain dynamical trajectories that are not stationary, 
thus obviating Einstein's concern, and allowing for the 
possibility of satisfying the correspondence principle. 

For purposes of this discussion, the approaches of 
Floyd and of FM may from this point forward be re- 
garded as identical. Since both employ the QSHJE, as 
does the present approach, it is clear that all three utilize 
essentially the same bipolar functions, s{x), r{x), p(x), 
and q{x). However, there are some important differences. 
Floyd uses a different convention, for which both the nor- 
malization of ^'(a;), and the flux, change simultaneously, 
so that his r{x) is proportional to ours. Neither of the 
other methods identifies the two parameters associated 
with Eq. (fT4|) in the way that we have done, and they 
certainly do not provide a means of selecting preferred 
values for these parameters, as per Sec. IIII C 31 

Indeed, Floyd considers each of the two-parameter 
family of solutions to Eq. ([H]), which he terms "mi- 
crostates," to constitute an equally valid decomposition 
of ^'(x). He further asserts that since the Schrodinger 
equation per se provides no means of distinguishing mi- 
crostates, that the QSHJE must be regarded as the more 
fundamental equation, in some sense. Floyd also ac- 
knowledges the fundamental differences between q{x) and 
and the advantages of the former vis-a-vis nodes, 
as discussed in Sec. Ill CI (although he incorrectly claims 
that the latter is singular at the origin for the first excited 
harmonic oscillator eigenstate^) . On the other hand, the 
lack of a criterion for selecting a prefered microstate im- 
plies that the quantum trajectories associated with his 
approach — termed "Floydian trajectories" — will not in 
general approach classical trajectories in the correspon- 
dence principle limits. 

Actually, there is another, more fundamental reason 
why Floydian trajectories do not satisfy the correspon- 
dence principle. This is because the underlying dynam- 
ical law governing their evolution is radically different 
from that of the Bohmian quantum trajectories utilized 
here. The trajectories are different, even though the bipo- 
lar decompositions are identical — a very curious situation 
that bears further analysis. In particular, the period of 
Floydian trajectories is finite, necessitating a speed that 
approaches infinity as — > oo. In contrast, Bohmian 
trajectories rapidly slow down in the forbidden regions, 
never reaching the coordinate asymptotes, as discussed in 
Sec. IIV Al Whereas Floyd finds implications for hidden 
variables and the Copenhagen interpretation, we adopt 
a decidedly less philosophically ambitious perspective. 

The Floydian dynamical law is obtained from the 
QSHJE by applying a classical procedure, wherein tra- 
jectory evolution is related to the quantity ds/dE to yield 



where u{x) = V{x) + q{x) is the bipolar modified po- 
tential. Various interpretations may be provided for the 
(l—du/dE) factor in Eq. ([26]) . FM adopt a relativistic in- 
terpretation which lumps it together with m to form the 
"effective quantum mass," whereas Brown (in essence) 
uses it to define an "effective time." In the Bohmian 
approach, u{x) is regarded as independent of E, thus 
giving rise to the usual p = mx relation. In the Floydian 
approach however, u{x) is considered to depend on E, 
giving rise to the more complicated Eq. (j26p expression 
above, for which the canonically conjugate momentum, 
p{x), is not necessarily equal to mechanical momentum, 
mi. 

How can the same u{x) be regarded as energy- 
independent in one theory and energy-dependent in an- 
other? It has to do with different interpretations of the 
energy. One can make a distinction between the quan- 
tum energy, Eq , and the classical energy, Ec ■ The quan- 
tum energy is that of the eigenstate ^E", which determines 
q{x), u{x), and h(x,p) — p^ / 2m + u{x) . The contours of 
h{x,p) define the classical energies Ec, with Eq — Eq 
corresponding to the bipolar quantum LM. The function 
u{x) clearly depends on Eq, but does not depend on Eq. 
The meaning of Eq. ((26| therefore depends on whether 
E = Eq (Floydian approach), or E = Eq (Bohmian ap- 
proach) . 

Floydian dynamics, therefore, represents a different 
philosophical outlook than that of Bohmian dynamics. 
Both approaches are correct, and ultimately yield math- 
ematically identical results. The main point we wish to 
make here is simply that this difference does not ap- 
pear to have anything to do with the bipolar ansatz per 
se. Indeed, one could easily apply the Floydian dynam- 
ical law to the unipolar ansatz — although for stationary 
eigenstates, this would yield the same results as Bohmian 
propagation. 

We close this section with a few final comparisons be- 
tween the Floyd/FM and present approaches. First, we 
comment that quantum energies are discrete for bound 
states, and so any derivative with respect to Eq requires 
careful consideration, as has been previously noted. ^^'^^ 
This is relevant for Floydian trajectories, but is essen- 
tially a non-issue for Bohmian trajectories. Second, al- 
though Floyd has certainly considered the 1 DOE har- 
monic oscillator system, he appears not to have derived 
analytic expressions for the relevant bipolar functions, as 
we have done in Sec. lIVI We have, moreover, verified that 
Floyd's closed form expression for the bipolar modified 
potential)^ 

u{x) =E-l/ [a-^i{xf + b^2{xf + c^i{x)^2{x)]^ , 

(27) 

is consistent with Eq. (|24p . by substituting 'Si = 
2rcos(s), 5*2 = 2rsin(s), a = b = c = 0, and 

E — 1/2 into Eq. ([27]) . Finally, we comment that 
the Floyd/FM approach, being based on the QSHJE, 
does not appear to generalize to the non-stationary-case, 
whereas the present approach does — at least in certain 
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situations. 



VI. SUMMARY AND CONCLUSIONS 

The Schrodinger equation is linear, yet the equivalent 
Eqs. ([4]) and ([5]) — obtained via substitution of the Eq. ([2|) 
ansatz into the Schrodinger equation — are not. Quite 
apart from the philosophically intriguing issues which 
this raises, a primary conclusion of the present work is 
that this situation may also be exploited for practical 
purposes. In particular, the superposition principle al- 
lows us to divide up the initial wavepacket "^(x) into 
pieces ^^(a;), evolve each of these separately over time, 
and then recombine them to construct the time-evolved 
^(x) itself. So far as the Schrodinger equation is con- 
cerned, the division into ^'^(a;) pieces is arbitrary. How- 
ever, if the evolution of the ^kix) is performed using 
Eqs. ([U and (O, nonlinearity implies that the division 
is not arbitrary, but has a large impact on the time- 
dependent behavior of the resultant Rk(x) and Sk{x). 

In principle, therefore, one may improve the numeri- 
cal performance of quantum trajectory calculations sim- 
ply by judiciously dividing up the initial wavepacket into 
several pieces. It remains to be seen the extent to which 
such a procedure will prove beneficial for actual numer- 
ical calculations of real molecular systems — to be sure, 
much depends on the manner in which '^(x) is decom- 
posed. Nevertheless, the particular bipolar scheme ex- 
plored in this paper, already appears to exhibit much 
promise — at least with regard to ameliorating the infa- 
mous node problem, which has thus far severely limited 
the effectiveness of QTMs in the molecular arena. 

The basic idea is to decompose a wavepacket, ^'(a;), 
that has nodes, into a linear combination of two com- 
ponents, 5'±(x), that do not. In practice, it is not only 
nodes per se that cause problems for numerical QTMs, 
but more generally, any large or rapid oscillations in 
^(x). Thus, ^±{x) should ideally be not only node- 
less, but also smooth and slowly varying. For the special 
case where "^{x) is a stationary Hamiltonian eigenstate, 
the semiclassical method is well-known to yield approx- 
imate ^±(x) functions with the requisite properties — 
even when ^'(x) itself is highly oscillatory. It is for this 
reason that the semiclassical solutions were used as a 
guide for determining the corresponding exact quantum 
^'-l-(x)'s. Not only have we provided an explicit recipe for 
obtaining the latter, we have also shown that these satisfy 
the correspondence principle in the appropriate semiclas- 
sical limits. Thus, the bipolar quantum potential, q{x), 
obtained here — now understood to represent the quan- 
tum correction to the semiclassical approximation — is 
not only well-behaved in the vicinity of nodes, but ac- 
tually approaches zero, in the large action limit. 

The bipolar quantum trajectories also behave very 
much like classical/semiclassical trajectories; indeed, the 
two are nearly identical in the classically allowed region, 
in the large action limit. This may seem paradoxical, as 



both bipolar and unipolar quantum trajectories conform 
to the same dynamical law, and the latter are known to 
behave very non-classically. Unipolar trajectories do not 
cross in position space, for instance — which can cause 
kinky trajectories and other node-related difhculties, es- 
pecially when the wavepacket undergoes reflection. The 
bipolar trajectories get around this difficulty as follows: 
whereas the trajectories on a single LM sheet never cross 
each other, they are all headed in the same direction 
anyway, and so they don't get in each other's way. On 
the other hand, trajectories on one LM sheet are free 
to cross those on the other sheet — ^just like the corre- 
sponding semiclassical trajectories. From a philosophi- 
cal standpoint, one might thus regard the present bipo- 
lar decomposition to be more compelling than the stan- 
dard unipolar approach — although curiously, this stance 
would require one to abandon Bohm's original pilot wave 
interpretation. 

The above discussion anticipates future application of 
the present ideas to arbitrary time-evolving wavepackets; 
but it must be borne in mind that thus far, only station- 
ary wavepackets have been considered. It is encourag- 
ing that the bipolar decomposition scheme outlined here 
was found to be preserved over time for stationary states 
(Sec. nil C ip . On the other hand, for the more general 
time-evolving case, the LMs themselves will change over 
time, and the bipolar decomposition scheme itself need 
not be preserved. It is therefore possible that the initially 
nodeless ^±(a;)'s may develop nodes over the course of 
time. This need not cause difhculties in practice however, 
because at any desired time, one is free to redecompose 
^(x) into new \E'±(x)'s that are nodeless. 

One of the most appealing aspects of the excited har- 
monic oscillator results of Sec. IIVBI is the fact that the 
bipolar functions remain smooth and slowly varying for 
all values of n. Indeed, the LMs for all n values re- 
semble each other, apart from a change of scale. This 
is very advantageous from a numerical perspective, as 
it suggests that very few trajectories would be needed 
to accurately compute r„(x) derivatives, if a completely 
numerical propagation scheme were adopted. More to 
the point: the number of trajectories required should be 
essentially independent of the number of nodes. For suf- 
ficiently large n, it should even be possible to perform 
an accurate calculation with fewer than n trajectories — a 
prospect that would be virtually unheard of in a unipolar 
context. Note that since the bipolar trajectories them- 
selves are also much smoother than the unipolar trajecto- 
ries, far fewer time steps should be required in the bipolar 
case. 

We conclude with a brief discussion of the prospects for 
multidimensional systems. At present, it is not entirely 
clear how best to apply bipolar decomposition to an arbi- 
trary, real, multidimensional wavepacket, ^(xi, . . .). For 
regular systems, semiclassical theory suggests an essen- 
tially direct-product decomposition, via pairs of action- 
angle coordinates. Generally speaking, however, it is dif- 
ficult to find these coordinates — unless the wavepacket is 
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initially separable, which is very often the case in molec- 
ular applications. On the other hand, the factor-of-two 
bifurcation applies to each degree of freedom separately, 
resulting in 2^ components total, where D is the num- 
ber of DOFs. This is clearly undesirable for large D. A 
more effective strategy may be to bifurcate '^{x) along 
the reaction coordinate only. These and other ideas will 
be explored in future publications. 
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